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ABSTRACT 

All estimation has been made of the principal long period spherical harmonic 
parameters in the representation for tlie M 2 ocean tide from the orbital histories of tlirec 
satellites — 1967-92A (TRANSIT), Starlettc, and GEOS-3. The data used were primarily 
the evolution of the orbital inclinations of the satellites, with the addition of the longitude 
of the ascending node from GEOS-3. The results are: 

C ^2 = 3.42 ± 0.24 cm. 42 

These values agree quite well with recent numericai models and another recent determina- 
tion from satellite data. 

Dissipational tidal friction in the oceans is known to provide the largest contribution to 
tlie observed deceleration in the lunar mean longitude. Further, only the second degree 
components of the ocean tide contribute significantly to this secular decay (Lambcck, 



1975). The M, pavnmeters obtained here infer an ^ of-25 ± 3 arc seconcls/ccntury^, in 
yood agreement with other investigators. The range of current determinations of is from 
-24.6 to -27.2 arc seecnd/century^. Considerably different techniques have been used to 
derive the estimates; the study of ancient eclipses, transits of Mercury, Lunar Laser Rang- 
ing, and another satellite solution. 


OCHAN TIDE PARAMETERS AND THE DECELERATION OF THE 
MOON’S MliAN LONCrrUDE FROM SA'rELLITE ORBH’ DATA 


I NTRODUCTION 

In the liist decade and a half, there have been a number of iiivcsligalions concorniiig the 
clTccts ofsolid and niiid earth tides on satellite orbits. With improved tracking tccliniques, 
more accurate modelling of the Earth’s gravity field, and more data from satellites in differ- 
ent orbiUs, it has been possible to extract increasingly better information from the satellite 
orbit data about tidal effects and the resulting iirplications regarding Umar orbit evolution. 

Ocean tidal effects are most easily obsen'ed in the inclination histories of satellites, it 
being the best behaved of the orbital elements for tidal studies, in some instances, the per- 
turbations can be seen satisfactorily in the longitude of the ascending node. Tidal constitu- 
ents such as Sa, Ivj , Ko, and P[ typically cause variations up to several tenths of an arc 
second in amplitude, witli periods of tens or hundreds of days. The most intriguing Ocean 
tide to study is the M^ because it has the greatest effect on the Earth’s oceans. However, Its 
influence on satellite orbital elements is small, usually no more than a few hundredths of an 
arc second in amplitude with periods averaging about two weeks. As a consequence, very 
precise data must bp .derived in order to obtain a significant recovery of parameters for the 
M 2 ocean tide. Table 1 presents a compaiison of the expected and obscivcd M 2 effects on 
the orbits of the three satellites msed in this investigation. 

Lambeck ct al. (1974) presented analytical cxprc.ssions in terms of spherical harmonics 
for the ocean tide effects on satellite orbits. They indicate that the obseivable perturbations 
arc long period in nature (on tlic order of days) and arc produced by only the first few low 
degree and order terms in the development of each tidal constituent. 

This formulation of Lambcck et al, (1974) has been adopted by a number of recent 
investigations of satellite derived ocean tidal effects. Fclsentroger ct al. (1976) and 



Folscnti-cgcrot al. ( 1978) estimatcii ocean lidc panimclers from GIiOS*l, G EOS-2, and 
1967-92A obscivational data, These were single satellite solutions making use ofindination 
data only, so .separation between the second and fourth degree harmonics could not be ob- 
tained. Cazenave ct al. (1977) studied the orbits of GEOS-1 and I967-92A and presented 
estimates of (he (2,2) and (4,2) parameters tor the So and Mo ocean tides. However, tiicir 
results for So arc probably corrupted by questionable solar radiation pressure modelling and 
their Mo values arc dominated by the 1967-92A data. Goad and Douglas (1977a) success- 
fully derived values of the (2,2) and (4,2) parameters for the M 2 ocean tide from 1967-92A 
and GEOS-.3 data. Tlicy also used these results to evaluate the deceleration of the Moon’s 
mean longitude, according to tlic procedure described by Lambcck (1975), 

In this i)resent work, we have used considerably more data for 1 967-92A and GEOS-3 
than either our earlier work or that used by Goad and Douglas (1977a),* Furtherm^e, to 
strengthen the solution we have incorporated data from Starlctte, wliich has a signirtbantly 
different inclination than the other two satellites. r 

ANALYSIS PROCEDURE AND ANALYTIC FORMULATION 

Our analysis procedure begins with obtaining a definitive set of satellite osculating 
orbital elements from the tracking data, These elements arc then converted to “mean” 
Kcplcrian orbital elements. The “meaning” procedure consists of removing the first order 
perturbations from the osculating elements, and then passing these preliminaiy mean ele- 
ments through a low-pass filter to numerically remove residual high frequency effects. The 
final mean elements are the result of llltcring over two day spans. Note that prior to (he 
nitcring, the direct lunar perturbations are removed in order to eliminate errors in the 
lunar ocean tide parameter determination caused by damping the direct effect due to the 

*Tiicrc has been a coasi(lcr.iblc amount of cooperation between Goad and Douglas and the authors for the 
purpose ofgetiorating (he precision mean element data. Thus, there is a substmitial data overlap, 
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averaging (Go!k 1 and OougUis. 1977a)» After nitcring, the efrccts of the direct lunar pertur- 
bations arc reinserted. 

Tlic mean elements arc then processed tluough an orbit determination program wliicli 
computes a theoretical orbit based upon long period, resonant and secular perturbations. 
This theoietical orbit is differenced with the mean eiements to produce a set of “residuals” 
from which the tidal perturbations can be isolated, 

In order to separate the ocean and solid Earth tidal effects, a solid Earth tide model 
with k 2 = 0.3 and lag angle 63 = 0° is used in the orbit determination procedure. This l <2 

I 

value has been well established to ±.01 from seismic data, Earth tide measurements, and 
Earth rotation obscivations. This modelling is necessary because the same frequencies arc 
present in both tides; otherwise, meaningful ocean tide parameters could not be separated. 
In addition, Goad and Douglas (1977b) have shown that the only satellite derived M 2 ocean 
tide parameters consistent with both the recent luunerical ocean tides and values of the de- 
celeration of the lunar mean longitude are those obtained with = 0.3 and 62 == 0 . 


The fonniilation we use for the ocean tide potential is expressed in terms of the equa- 
torial Kcplcrian orbital elements of the satellite (Lambcck ct al., 1974; Felscntregcr ct ah, 
1976) 

P j=0 p=0 <!=-« ^ 


where 




j- fi oven 
j-fi odd 


(I) 


T'fcpq " ^ ^ " ^9 + 9)M + “ 0) i 27rfr ± 
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g = menu acceleration of gravity at surface of Earth 
p = mean tlcnsity of sea water 
p = average density of the Earth 
kj = ioad deformation coefficient of degree j 
R = mean radius of the Eartii 
a = semi-major axis 
i = inclination 
c = eccentricity 
= argument of perigee 
M = mean anomaly 
= longitude of ascending node 
0 = Greenwich mean sidereal time 

“ inclination and eccentricity functions (Kaula, \9C6) 

2irCT = argument for tidai constituent 
Cjj, = spherical harmonic coefficients 
erg = phase angles 

Long period terms occur only with Cjg and when j - 2p + q = 0, Perturbations in the 
elements of each satellite can be found by substituting (1) into the Lagrange Planetary 
Equations and performing the integration under the assumption that only the arguments are 
time varying. The resulting perturbation equations for the principal second and fourth de- 
gree harmonics of the global Mo ocean tide are as follows: 

1967-92A 


A r _ 0.0 1 003 ai C sec /n r-\ aia / r _+ A 

At Co 2 sm (2f2 - 2\ + £ 22 ) 


• 22 ^ 


+ 


0.009606 arc sec . 


cm 


CJ 2 sin (2f2 - 2\' + £^ 2 ) 


( 2 ) 
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STARLIi'H'E 


GEOS-3 


AI = 


0.006268 arc sec ^ . . + 

— CI 2 sin (2n - 2X + £ 22 ) 


cm 


0.01 190 arc sec ^ 4 . . ^ 

■ — CJ 2 sin ( 2 S 2 - 2 \' + 


cm 


42 ^ 


0.01288 arc sec ^ 4 . . ^ 

C+ 5IH - 2X + et^j) 

cm 


0.0032j 7 aro je c * ^ , 

cm ^ 


( 3 ) 


(4) 


-0,002394 arc sec _j. ... + . 

AS2 Cto cos (2fi - 2X^ + 

cm 


0,03464 arc sec 
cm 


CJ 2 cos (2i2 - 2 V + 42 ) 


A2> 


(5) 


where 

X' = mean longitude of the Moon (referred to ecliptic) 

At this point, we correct for the two-clay averaging, When the averaging was performed, 
both the solid Earth and ocean tidal effects were present in the data. Therefore, in order to 
recover tlie uncomipted ocean tide signals, the solid Barth tide model used in the orbit de- 
termination step must be reconstructed and a correction applied. The solid Earth tide 
potential outside the Earth is (Lambeck, 1975) 


C=2 111=0 [)=0il=-“ j=0 ^ / \ / 


^*^Cnipq ~ '^Sinjli ‘^.Cinpq^ 


( 6 ) 


where 
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= (fi~ 2p)o>* H* (fi- 2p + q)M* + mfi* 

- 2j)co + (fi - 2j + lOM + mS2 

Here, (ii*, e*, I*, M*, to*, fi*^) arc llic equatorial Kcplerlan elements of the Moon, in* is 
the mass of the Moon, and kg Is the Love number of degree B. 

Also, 




1 , m=0 

2, niT^O. 

Again, the long period part of the potential is extracted and substituted into the Lagrange 
Planetary Equations which are integrated as before. The perturbation equations which arc 
generated are the following (for the M^ constituent): 

1967-92A 

AI[j = (1,04713 arc sec) k2 cos (*^2200 " *^2200^ 


ST ARLETre 

AI|^ “ (0,6dS13 arc sec) ko cos (^^7200 ^ 2f2 4* ^'2200^ 

GHOS-3 

AIg = (1.36861 arcsec)ko cos(r'2200 ~ 2fi + e2200^ 
Af2j7 = (-0.63760 arc sec) k2 sin (t-Jjoqo “ 20 + ^2200^ 


( 8 ) 

(9) 

( 10 ) 


OCEAN tidal parameter ESTIMATION 

Data utilized for the 1967-92A, Starlcttc, and GEOS-3 satellites spanned time periods 
of 3 12, 252, and 382 days, respectively. Frequency analyses were performed on the residuals 
to ascertain the periodicities to be included In a series of least squares fits to the residuals. 
These least squares analyses produced amplitudes and phases for the observed harmonic vari- 
ations in the residuals. Figures 1 , 2, 3, and 4 display the residuals and the fits, along with 
the observed amplitudes for the M-, effects. The oscillations caused by the M2 tide can be 
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clearly seen in each instance. The mis of fits for the four figures arc .036, .051, .030, and 
,034 arc seconds, respectively. 

The expressions in Equations 2, 3, 4, 5, 7, 8, 9, 10, along with the ainpiitudes and 
phases for the Mt title obtained in the breakdown of the residuals, were then used in a sub- 
sequent least squares solution for tlie (2,2)'*' and (4,2)^ Mo ocean tide parameters, In the 
process, the correction for the averaging was made (see Jeffrcy.s aiio Jeffreys, 1972) and tlic 
values k^ = 0.3 and 82 = 0“ were used for the solid Eartli tide parameters. The results are 
shown in Table 2, along with several other recent determinations. The uncertainties arc the 
formal standard deviations of the solution the realistic errors are probably somewhat 
larger. 

DECELERATION OF MOON*S M E AN LONGITUDE 

The tidal potentials (1 ) and (6) can be used to study the effects of the tides on the 
lunar orbit evolution. Lam beck (1975) pointed out that the tidal perturbation expressions 
obtainable from these potentials arc valid for any satellite of the Earth, and thus arc applica- 
ble to the Moon itself. Me also demonstrated that the deceleration of the lunar mean longi- 
tude was caused principally by the second degree components of the ocean tide and possibly 
the second degree solid tide (if the phase angles arc assumed to be different from zero). A 
more recent paper by Lambeck (1977) gives an c.xccllent historical and technical oveiview 
on the subject of the dissipation of tidal energy in the oceans. 

The mean motion N of the Moon can be expressed as: 



where M is the mass of the Earth. Taking the time derivative of (! 1) results in 
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wlileli gives the secular rate of change of N (Lambcck, 1975; 1977). The purely secular 
tenns in (6) arc those which “ ''cnijli " or p = j and q = h. These tcnns arc sub- 
stituted into the Lagrange Planetary Equation for a, resulting in an expression fora which 
is inserted in (12), Then, the contribution to N from the solid Earth tide is 


8 fi « 


nb- z 11 l“i ra 

lia3 ninfl non nn-M ' / 


l^\2C+l(g- 


(fi + ni)l 
* (^ “ 2p + q) sin Cgpprt 


8“2 m°0 i»=D <!= 
|2 


(13) 


Table 3 suinniarizes the principal second degree solid tide contributions to N. 

The main ocean tide variations in N are provided by the second degree tcnns for the 
M2, N2. and Oj constituents, and only for Tj^pq (Lambcck, 1 975). In the ocean tide 
potential (I), the (nearly) secular M., term arises when j = 2, J2 = 2, p « 0, q = 0; for N2, j = 2, 
J2 = 2, p = 0, q = 1; for 0| , j " 2, 1, p = 0, q = 0. Tabled presents these contributions. 

In order to compute a value for N, we use the fact that the solid Earth tide phase 
angles are approximately equal to nifio, where ^2 is the solid Earth tide lag angle. In 

this work, the assumption has been made that 5^ = 0°. Goad and Douglas (1977b) have 
sirown that a change in the assumed values for k2 and 62 will result in a corresponding alter- 
ation in the ocean tide parameters derived from satellite data. Tiicrcfore, tlie solid Earth 
tide contributions to N given in Table 3 must be taken to be zero. 

This leaves the ocean tidal effects listed in Table 4. \Vc use our satellite derived M2 
parameters presented in Table 2 and the N2 and O j models of R, Estes of Business and 
Technological Systems, Inc,, Senbrook, Maryland: 


^2 



.58 cm : 358® 


O] 


'-21 

^21 

2.70 Cm 

159° 
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'rhe result is a value tor N of -25 ± 3 arc scconclF/ccntury^, Table 5 presents a comparison 
of this estimate with several other recent determinations, 

DiSCUSSiON 

The satellite derived (2, 2)'^ parameters for the M 2 ocean tide arc in accord with the 
numerical models of Estes and Schwiderski and the satellite solution of Goad and Douglas 
(1977a), if one assumes the realistic uncertainties to be at least 50% larger than the formal 
standard deviations. The (4,2)'*‘ values are in good agreement, with the noticeable excep- 
tions of the C 42 of Schwiderski and the 642 of Estes. The data used by Goatl and Douglas 
(1977a) form a subset of the data used in the present investigation, so the consistency be- 
tween the two models is not unexpected, ifowever, the agreement with the numerical 
models is quite satisfying. 

The credibility of oiir M 2 ocean tide solution is further enhanced by the close accord 
betw ^u v.itr computed value for the cicccicraflon of the lunar mean longitude and the other 
recently reported estimates. Under the assumption that k 2 = 0.3 and §2 “ O'* for the solid 
Eartli tide, the major contribution to the N comes from the (2,2)"^ parameters for the M^ 
ocean tide. 

It is evident from the results of this Investigation that :jtiidies of close Earth satellite 
orbits arc able to provide important information about the tidal forces acting on the Earth. 
We anticipate that the new laser data for Starlctte and GEOS-3, in addition to data for the 
SEASAT Satellite, will be even more preci.se for studies of this sort. 
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TAI3LI:3. SUMMARY OF PRINCIPAL SOLID EARTH TIDAL EFFECTS ON THE 
DECELERATION OF THE MOON’S MEAN LONGITUDE 


TIDAL CONSTITUENT 

AN (ARC SECONDS/CENTURY^) 

M 2 

-870.8 k 2 sill ^2200 

N 2 

““48«3 k sin ^220J 

Ol 

- 1 12,1 IC 2 sin C2100 

L 2 

-0.3 IC 2 sin e 220 (-i) 

•^1 

0.4 IC 2 sill e 2 j,(_i) 

Qi 

-6.2 IC 2 sin «? 2 ]oi 
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TABLE 4. SUMMARY OF PRINCIPAL OCEAN TIDAL EFFECTS ON THE 
DECELERATION OF THE MOON’S MEAN LONGITUDE 


TIDAL CONSTITUENT 

AN (ARC SECONDS/CENTURY“) 

M 2 

CI 2 (-8.04 cos cJt + 0-22 sin 

N 2 * 

C 20 (-2.32 cos £22 *^22^ 

Oi* 

C 21 (-0.32 coset 1 - 1.41 sin ejj) 

L_ ’ 

*Estes (BTS) numerical models (obtained from C. Goad, private communication). 
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TABLE 5. RECENT VALUES FOR THE DECELERATrON 
OF THE MOON’S MEAN LONGITUDE 


SOURCE 

N (ARC SECONDS/CENTURY^) 

FELSENTREGERctaL* 

“25 ± 3 

GOAD AND DOUGLAS ( 1 911h)- 

“26 ± 3 

MORRISON AND WARD (1975)3 

“26 ±2 

MULLER (1976)'^ 

-27,2 ±1.7 

CALAME AND MULHOLLAND (1978)5 

-24.6 ±1.6 


* Satellite Solution, 

^Satellite Solution, 

^From Transits of Mercury. 

‘^Froin Ancient Eclipses and Transits of Mercury. 
^Froin Lunar Laser Ranging. 
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Figure 1. Inclination Residuals for 1967^92A (Transit) 


M2 period = 10.5 DAYS 



Figure 2. Residuals in Inclination for Starlette 
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igure 3. Inclination Residuals forGEOS-3 



